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1.  INTRODUCTION: 


In  medical  studies  or  in  reliability  analysis  it  is  quite  common  that  more  than  one  risk 
factor  may  be  present  at  the  same  time.  An  investigator  is  often  interested  in  the  assesment  of 
a  specific  risk  in  presence  of  other  risk  factors.  In  the  Statistical  literature  it  is  well  known  as 
the  analysis  of  competing  risks  model.  In  analyzing  the  competing  risks  model  it  is  assumed 
that  the  data  consists  of  a  failure  time  and  an  indicator  denoting  the  cause  of  failure.  Several 
studies  have  been  carried  out  under  this  assumption  for  both  parametric  and  nonparametric 
set  up.  For  the  parametric  set  up  it  is  assumed  that  the  different  lifetime  distributions 
follow  some  special  parametric  distribution,  namely  exponential,  gamma  or  Weibull.  Several 
authors  for  example  Berkson  and  Elveback  (1960),  Cox  (1959),  David  and  Moeschberger 
(1978)  considered  this  problem  from  the  parametric  point  of  view.  In  the  nonparametric  set 
up  no  specific  lifetime  distribution  is  assumed.  Kaplan  and  Meier  (1958),  Efron  (1967)  and 
Peterson  (1977)  analysed  the  nonparametric  version  of  this  model.  But  in  all  the  above  cases 
it  is  assumed  that  the  causes  of  failure  are  known  if  the  failure  times  are  observed.  But  in 
certain  situations  (Dinse;  1982  or  Miyawaka;  1982)  it  is  observed  that  the  determination  of 
the  cause  of  failure  may  be  expensive  or  may  be  very  difficult  to  observe.  Therefore  it  might 
occur  that  the  failure  time  of  that  item/  individual  is  observed  but  the  corresponding  cause  of 
failure  is  not  observed.  Miyawaka  (1984)  considered  this  model  and  obtained  the  maximum 
likelihood  estimators  (MLE’s)  and  the  uniformly  minimum  variance  unbiased  estimators 
(UMVUE’s)  of  the  failure  rates  of  the  different  failure  distribution  under  the  assumption 
that  they  are  of  the  exponential  type.  But  he  did  not  provide  any  distributional  properties 
of  these  estimators. 

In  this  paper  we  consider  the  same  model  as  that  of  Miyawaka  (1984).  It  is  assumed  that 
every  member  of  a  certain  target  population  either  dies  of  a  particular  cause,  say  cancer,  or 
by  other  causes.  A  proportion  tt  of  the  population  die  of  cancer  and  the  proportion  (1  —  7r) 
die  due  to  other  causes.  At  the  end  of  the  study  ,  we  have  three  types  of  observations: 


(a)  Individuals  who  died  of  cancer  and  their  lifetimes. 

(b)  Individuals  who  died  of  other  causes  and  their  lifetimes. 

(c)  Individuals  who  died  of  unknown  causes  their  lifetimes  are  observed,  but  their  causes 
of  death  are  not  known. 

There  is  another  dimension  to  this  problem.  The  research  project  is  financed  for  a  fixed 
length  of  time,  say,  M.  Some  individuals  could  be  alive  at  the  end  of  the  project  period.  For 
simplicity,  we  ignore  this  aspect  of  the  problem  and  assume  that  every  one  in  the  sample  can 
be  monitored  until  theirs  death.  We  assume  that  the  lifetime  distributions  of  the  different 
causes  of  failure  are  exponential.  It  is  observed  that  although  the  MLE’s  or  the  UMVUE’s 
of  the  hazard  rates  always  exist  but  the  MLE’s  or  the  UMVUE’s  of  the  mean  lifetime  of 
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the  different  causes  may  not  exist  always.  It  is  not  possible  to  obtain  the  UMVUE’s  of 
the  mean  lifetimes  of  the  different  causes.  We  obtain  the  conditional  MLE’s  of  the  mean 
lifetime  of  the  different  causes.  We  obtain  the  exact  distribution  of  the  conditional  MLE’s 
using  the  conditional  moment  generating  function  approach  and  also  obtain  the  asymptotic 
distributions  of  the  MLE’s.  It  is  observed  that  the  MLE  of  the  mean  lifetime  usually  over 
estimate  the  true  parameter  at  least  for  small  sample.  As  in  complete  sample,  the  bias  can 
be  obtained  as  the  inverse  moment  of  the  positive  binomial  random  variable  even  in  the  case 
of  incomplete  data.  We  provide  small  tables  for  the  biases  of  the  MLE’s  for  different  sample 
sizes  and  for  different  parameter  values  when  10%  and  20%  data  are  incomplete.  Based  on 
the  exact  and  the  asymptotic  distributions,  we  propose  two  approximate  confidence  intervals 
of  the  different  parameters  of  interest.  We  also  use  percentile  bootstrap  and  bootstrap-t 
confidence  intervals  for  the  unknown  parameters  and  compare  their  performances  through 
Monte  Carlo  simulatons. 

Since  the  exponential  distribution  has  constant  failure  rates  it  might  not  be  very  practical 
to  assume  that  the  lifetime  distribution  is  exponential.  Two  parameter  Weibull  distribution 
can  be  used  to  analyze  the  lifetime  data  because  of  its  increasing  and  decreasing  failure  rates. 
We  consider  the  model  when  the  underlying  lifetime  distributions  are  Weibull.  We  obtain  the 
maximum  likelihood  estimators  of  the  different  parameters  and  study  their  properties  under 
this  assumption.  We  also  obtain  the  asymptotic  confidence  bounds  and  the  bootstrap  con¬ 
fidence  bounds  of  the  different  parameters  and  compare  their  performances  through  Monte 
Carlo  simulations.  We  consider  one  data  set  from  Lawless  (1982)  and  see  how  the  different 
methods  work  in  this  practical  situation. 

The  rest  of  the  paper  is  organised  as  follows.  In  Section  2,  we  give  the  description  of  the 
model  and  give  the  different  notations  we  are  going  to  use  in  this  paper.  In  Section  3,  we 
consider  the  estimation  of  the  different  parameters  and  also  obtain  the  exact  distribution 
of  the  MLE  of  the  mean  lifetime  when  the  lifetime  distribution  is  exponential.  Different 
confidence  intervals  of  the  unknown  parameters  are  considered  in  Section  4.  Weibull  lifetime 
distribution  is  considered  in  Section  5.  The  numerical  results  are  presented  in  Section  6.  One 
data  set  from  Lawless  (1982)  is  being  analysed  in  Section  7  and  finally  we  draw  conclusions 
from  our  work  in  section  8. 


2.  MODEL  DESCRIPTION  AND  NOTATIONS: 

Without  loss  of  generality  we  assume  that  there  are  only  two  causes  of  failure.  We  use 
the  following  notations: 

X{\  Lifetime  of  system  i. 

Xji :  Lifetime  of  mode  j  of  system  i,  j  =  1,2. 
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F(.):  Cumulative  distribution  function  of  X{. 

Fj(.):  Cumulative  distribution  function  of  Xji,  j  =  1,2. 

m  =  i  -  m 

5ji  indiactor  variable  denoting  the  cause  of  failure  of  system  i. 

I[.]:  indiactor  function  of  event  [.]. 

Gamma  (a,  A):  denotes  the  gamma  random  variable  with  density  function  pja:Ql“1e_Al 

Weibull  (a,  A):  denotes  the  Weibull  random  variable  with  density  function  ot\xa~le~Xx 

It  is  assumed  that  (Xij,X2j);  i  =  1,2,  n,  are  n  independent  identically  distributed 
(i.i.d.)  random  variables.  Ah  and  Xu  are  independent  for  all  i  =  1,2,  . . .,  n  and  Xi  =  min 
{Xh,X2i}.  With  out  loss  of  generality  it  is  assumed  that  the  first  m  observations  consist  of 
failure  times  and  also  causes  of  failure  whereas  for  the  last  (n-m)  observations  we  only  observe 
the  failure  times  and  not  the  causes  of  failure,  i.e.  the  following  data  are  being  observed 
{Xu  <50,  (X2, 52), . . . ,  (Xm,  5m),  {Xm+U  {Xnt  *)•  In  order  to  analyze  the  incomplete 

data  it  is  assumed  that  the  failure  times  are  from  the  same  population  as  the  complete  data, 
that  is  population  remains  unchanged  irrespective  of  the  cause  of  failure. 


3.  EXPONENTIAL  FAILURE  DISTRIBUTIONS,  ESTIMATION: 


In  this  section  we  assume  that  XjiS  are  exponential  random  variables  with  parameters 
A  j  for  i  =  1,  2,  . . .,  n  and  for  j  =  1  and  2.  The  distribution  function  Fj{.)  of  Xji  has  the 
following  form: 

Fj(t)  =  (1  —  (3.1) 

for  j  =  1  and  2.  The  likelihood  function  of  the  observed  data  (xi, <Ji), . .  •  (xm,  5m),  (xm+i,  *), 
. . .,  (xn,  *)  for  the  general  case  takes  the  following  form: 


m  n 

L=n[c/F1(xi)F2(^)]/(<5i=1)[dF2(xi)F1(xi)]^=2)  dF(xJ-  ^ 

i=l  i=m+l 

Therefore  for  the  particular  case  when  F\  and  F2  are  exponentials  with  parameters  Ai  and 
A2  respectively,  the  likelihood  function  becomes: 


L  =  Ai'A^A!  +  A2)n  mexp{-{\i  +  X2)^2xi).  (3.3) 

t=i 
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Here  r\  and  r2  denote 
the  logarithm  of  (3.3) 
Aj  as 


the  number  of  failures  due  to  mode  1  and  mode  2  respectively.  Taking 
and  equating  the  partial  derivatives  to  be  zeros,  we  get  the  MLE’s  of 


A 


l  = 


nrx 


ra£"=iZ»’ 


(3.3) 


and  the  UMVUE’s  of  Ai  as 


\*  -  (n~  Du 
1 


(3.4) 


see  Miyawaka  (1984).  Therefore  MLE  of  the  survival  function  due  to  cause  1  (say  cancer)  is 


Fj{x)  =  e 


_ 


(3.5) 


and  the  MLE  of  the  hazard  rate  or  the  instantaneous  death  rate  due  to  cause  1  is  given  by 


dFi(x) 


—  Ai. 


F  i(x) 

The  realtive  risk  rate,  7r,  due  to  cause  1  is 

P[ Xu  <  X2i]  =  [°°  A Ie-X'xe-X'xdx  =  —  , 

and  because  of  the  invariance  property  of  the  MLE,  the  relative  risk  rate  7r  due  to  cause  1  is 


Ai 


7 r  = 


Ai 


Ai  +  A2 


n 

m 


(3.6) 


For  the  exponential  lifetime  distribution  as  (3.1),  Ai  represents  the  hazard  rate  and  0i  = 
denotes  the  mean  lifetime  due  to  cause  1.  Although  the  MLE’s  and  the  UMVUE’s  of  Ai 
always  exist  but  the  same  is  not  true  for  9i .  The  UMVUE  of  0\  does  not  exist  and  the  MLE 
of  6{  exist  only  when  rq  >  0.  The  conditional  MLE  of  #i,  say  9\,  when  n  >  0  is  as  follows: 


0i  = 


m  £?=i  Xi 


nr  i 


(3.7) 


If  n  =  0,  0i  does  not  exist.  Now  we  obtain  the  conditional  distribution  of  the  MLE  of 
9\,  conditioning  on  ri  >  0.  Our  approach  to  produce  the  confidence  bound  for  0i  is  based 
on  the  distribution  of  0i  similarly  as  Chen  and  Bhattacharya  (1988).  Moreover,  use  of  the 
MLE  provides  a  safeguard  against  any  obvious  loss  of  informations  and  ensures  asymptotic 
optimality  of  the  present  method.  Consider  the  following  lemmas: 

Lemma  1:  The  conditional  moment  generating  function  (mgf)  of  0i,  ^(t),  is  of  the  fol¬ 
lowing  form 

*!,(«)  =  £[e“'|ri  >0] 


5 


m\  - 1 


=  (1-0 


=  Eft  1- 


ml 


02 


0i 


^  i!(m  -  t)!  \0i  +  02)  \0\+02t 
tm0 1 0 


1  - 


tm0\02 
ni(0 1  +  02)  / 


l”2 


i=l 


ni{0 1  +  02) 


here 


q  = 


+  02 

for  i  =  1,  2  ,  . . m. 

Proof  of  Lemma  1: 


a„dPi  =  (l-0-‘  — 


$2 


i!(m  —  i)!  \0\  +  02 )  \9\  +  02, 


Note  that  £"=l  X{  is  Gamma  (n,  Ai  +  A2)  random  variable  and  r\  is  a  binomial  random 
variable  with  parameters  m  and  A^Aa . 

h,(*)  =  £[^'|ri>0] 

=  [e“‘  In  =  i]  P(n  =  i|rL  >  0). 

1=1  j 

Using  the  facts  p{  =  P(rt  =  i\n  >  0)  for  i  =  1,2,  . . .,  m  and  the  moment  generating  function 
of  Gamma  (a,  A)  is  (1  —  ^)~°\  the  result  follows  immediately. 

Theorem  1:  The  conditional  probabilty  density  function  (pdf)  of  0 1,  say  /^(x),  or  the 
conditional  probabilty  distribution  function,  say  becomes 

m  m 

/*, (x)  =  E Pi9i (x),  Fdi(x)  =  Y,PiGi{x), 

1=1  2=1 

where  gi(x)  and  Gi(x)  denote  the  density  function  and  the  distribution  function  respectively 
of  a  gamma  random  variable  with  shape  parameter  n  and  scale  parameter  f°r  i  = 

1,  2,  . . .  m. 

Proof  of  Theorem  1:  Obvious  from  Lemma  1. 

Therefore  the  conditional  distribution  of  0\  is  a  mixture  of  m  gamma  random  variables. 
It  may  be  also  noted  that  when  m  =  n,  we  get  the  distribution  of  0\  in  the  competing  risk 
model  when  the  lifetime  distributions  of  the  different  causes  are  exponential  and  there  is  no 
censoring  (see  David  and  Moeschberger;  1978).  With  the  help  of  Theorem  1,  we  can  obtain 
different  conditional  moments  of  0\.  We  give  the  first  and  the  second  conditional  momemts 
of  0\.  For  brevity  we  denote  them  as  E^)  and  E(0i)  respectively. 


E(0 1) 


m6 10  2 

(0i  +  02) 
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and 


E{6\) 


m29\9\n(n  +  1 )  ^  Pi 
{9\+92)2ti2  fr[i2' 


Since  in  both  the  cases  the  quantities  within  the  summation  sign  denote  the  inverse  moments 
of  the  positive  binomial  random  variable  it  is  not  possible  to  give  the  exact  expressions. 
Therefore  it  is  difficult  to  obatin  the  exact  bias  from  there.  The  tabulated  values  of  the 
first  moment  of  the  inverse  positive  binomial  random  variable  are  available  (see  Edwin  and 
Savage;  1954),  it  may  be  used  to  make  some  bias  correction.  It  is  not  pursued  here.  If  Z 
is  a  binomial  random  variable  with  parameters  N  and  P,  then  for  large  N,  E(^| Z  >  0)  « 
—  lh  and  E(^|Z  >  0)  «  (E(zj)*  =  (NP)7  see  Mendenha11  and  Lehmann  (1960).  Using 

those  approximatons,  we  can  say  that  for  large  m  and  n,  E(9X)  ~  9X  and  E{9\)  ~  9\.  It 
shows  asymptotically  9\  is  an  unbiased  estimator  of  9\ .  As  the  variance  of  9X  goes  to  zero, 
9X  is  a  consistent  estimator  of  9 1  also. 


We  present  two  small  tables  which  give  the  numerical  values  of  the  biases  for  different 
sample  sizes  and  for  different  vaules  of  92.  In  all  these  calculations  we  have  kept  9\  =  1. 
Table  3.1  represents  the  negative  value  of  the  biases  when  10%  of  the  data  are  incomplete  and 
Table  3.2  represents  the  negative  value  of  the  biases  when  20%  of  the  data  are  incomplete. 


Table  3.1 

Negative  biases  of  9{  when  10%  of  the  date  are  incomplete. 


n 

92  =  1.25 

92  =  1.50 

92  =  1.75 

92  =  2.00 

92  =  2.25 

92  =  2.50 

10 

.1284 

.1041 

.0868 

.0741 

.0646 

.0572 

20 

.0528 

.0431 

.0363 

.0315 

.0277 

.0248 

30 

.0330 

.0271 

.0231 

.0201 

.0177 

.0159 

40 

.0241 

.0198 

.0169 

.0147 

.0130 

.0117 

50 

.0189 

.0156 

.0133 

.0116 

.0103 

.0093 

Table  3.2 

Negative  biases  of  9X  when  20%  of  the  date  are  incomplete. 


n 

92  =  1.25 

92  =  1.50 

92  =  1.75 

92  =  2.00 

92  =  2.25 

92  =  2.50 

10 

.1469 

.1207 

.1011 

.0867 

.0756 

.0669 

20 

.0610 

.0495 

.0417 

.0360 

.0317 

.0284 

30 

.0377 

.0309 

.0263 

.0228 

.0202 

.0181 

40 

.0273 

.0225 

.0192 

.0167 

.0148 

.0133 

50 

.0214 

.0177 

.0151 

.0132 

.0117 

.0105 

From  the  Table  3.1  and  Table  3.2  it  is  clear  that  the  bias  of  the  MLE  of  9\  is  always 
negative,  that  means  9X  always  over  estimates  9\.  As  the  sample  size  increases  the  bias  also 
decreases.  If  we  have  more  incomplete  data  the  bias  is  also  more.  For  example,  if  9X  =  50, 
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and  =  1.50  and  only  90  %  of  the  data  are  complete  for  a  sample  of  size  40,  then  the  bias 
of  9 1  can  be  obtained  from  the  Table  3.1  as  40  x  .0198  =  .792.  Therefore,  it  is  not  a  very 
serious  bias. 


4.  EXPONENTIAL  FAILURE  DISTRIBUTIONS,  CONFIDENCE 
INTERVALS: 


In  this  section  we  propose  four  different  confidence  intervals  of  9\.  The  first  one  is 
based  on  the  distribution  of  9\  under  the  similar  kind  of  assumptions  as  that  of  Chen  and 
Bhattacharya  (1988).  The  second  one  is  based  on  the  asymptotic  distribution  of  9\.  We 
propose  to  use  percentile  bootstrap  confidence  interval  and  bootstrap-t  confidence  interval 
and  give  their  implementation  procedures  in  this  section. 

4.1:  Approximate  Confidence  Bound: 

First  let’s  assume  that  62  is  known.  Suppose  Pex[9\  >  b ]  is  a  monotonically  increasing 
function  of  9i,  and  let  b(0)  be  a  function  such  that  |  =  Pgx{9\  >  6(#i)]-  Then  for  9\  <  9[, 
we  have 

|  =  Pe[[9 1  >  b(9[)}  =  P$l[9i  >  b(9 1)]  <  P^i  >  b(9 0]. 

This  implies  that  b(9i)  <  b(9[),  that  is  b(6)  is  an  increasing  function.  Thus  b~l(9)  exists  and 
is  an  increasing  function.  So  we  get  1  -  |  =  <  9\  ],  it  implies  &_1(#i)  *s  the  lower 

bound  of  the  (1  —  a)  100%  confidence  bound  of  au.  If  9obs  denotes  the  observed  value  of  9\ 
find  9L  =  b~l(9obs)  such  that  §  =  P9l{9  1  >  9obs),  which  is  equivalent  to  find  9L  such  that  1 
-  f  =  PoL  (0t  <  9obs).  As  P9l[9i  >  6]  is  a  monotonically  increasing  function  of  #i,  therefore 
Pei[9i  <  c]  is  a  monotonically  decreasing  function  of  0\.  Let  c (9)  be  a  function  such  that  | 
=  PeAO  1  <  c(0i)].  It  can  be  shown  exactly  as  before  that  c (9)  is  a  decreasing  function  and 
therefore  c~l{9 )  exists.  Find  9u  =  c~l(9obs),  such  that  |  =  P9u(9i  <  9obs).  Since  the  closed 
form  expression  of  the  function  b(9)  or  c(9)  is  not  possible,  we  need  to  use  some  iterative 
technique  to  get  9i  and  9y.  Note  that  we  can  get  an  exact  (1  —  a:)  100%  confidence  bound 
on  6\,  if  we  know  92 .  Since  92  is  usually  unknown,  we  need  to  estimate  92  and  we  get  an 
approximate  (1  —  a)  100%  confidence  bound  on  0\.  The  construction  of  the  confidence  bound 
of  9y  is  based  on  the  assumption  that  P9l  [^1  >  6]  is  a  monotonically  increasing  function  of 
9.  It  is  quite  difficult  to  prove  this  assumption  because  of  the  complicated  nature  of  the 
function.  Some  heuristic  justification  can  be  given  as  follows.  Since  ^1  is  the  MLE  of  the 
mean  life  9X  due  to  cause  1,  for  fixed  02,  the  larger  the  parameter  9\  is,  the  more  probable 
it  will  be  for  its  MLE  to  exceed  a  given  value.  Numerical  values  of  Pox[9i  >  b]  for  various 
values  of  9\  and  b  confirms  this  conjecture. 
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Table  4.1 

The  tabulated  values  of  Pgx[9 \  >  6],  for  02  =  2,  b  =  1  and  when  10%  of  the  data  are 

incomplete. 


n  0i  = 

1.00 

1.25 

1.50 

1.75 

2.00 

2.25 

2.50 

2.75 

3.00 

10 

.475 

.685 

.815' 

.890 

.933 

.958 

.973 

.982 

.987 

20 

.483 

.766 

.905 

.962 

.985 

.994 

.997 

.999 

.999 

30 

.486 

.818 

.948 

.986 

.996 

.999 

1.00 

1.00 

1.00 

40 

.488 

.855 

.971 

.995 

.999 

1.00 

1.00 

1.00 

1.00 

50 

.489 

.883 

.983 

.998 

1.00 

1.00 

1.00 

1.00 

1.00 

Table  4.2 

The  tabulated  values  of  Pgx[9\  >  b ],  for  92  =  2,  b  =  2  and  when  10%  of  the  data  are 

incomplete. 


n  0i  = 

1.00 

1.25 

1.50 

1.75 

2.00 

2.25 

2.50 

2.75 

3.00 

10 

.045 

.130 

.246 

.371 

.486 

.585 

.667 

.732 

.784 

20 

.009 

.058 

.172 

.329 

.491 

.630 

.738 

.817 

.873 

30 

.002 

.028 

.125 

.297 

.492 

.661 

.785 

.868 

.920 

40 

.000 

.014 

.093 

.271 

.493 

.686 

.821 

.903 

.949 

50 

.000 

.007 

.070 

.249 

.494 

.708 

.849 

.927 

.966 

Tables  4.1  and  4.2  indicate  that  Pgx[9\  >  b ]  is  an  increasing  function  of  9\ 


4.2:  Asymptotic  Confidence  Bound: 

In  this  subsection  first  we  obtain  the  asymptotic  distribution  of  Ai  and  A2,  using  that 
we  obtain  the  asymtotic  distributon  of  9\  and  92.  The  Fisher  Information  matrix  of  the 
parameters  Xi  and  A2  is  I(Ai,  A2)  =  /»j(Ai,A2),  for  i  =  1  and  2.  Here 


and 


hj{  Ai,  A2)  — 


(dHogLj  At,A2)\ 
V  dXidXj  ) 


/ii(A,,A2) 


■Tl2(Ai,  A2) 
-^22 (Ai ,  A2) 


mX2  +  nX\ 
Ai(Ai  +  A2)2’ 

mAi  +  nA2 
A2(Ai  +  A2)2 


The  Fisher  Information  matrix  of  9y  and  92,  I(9i,92),  can  be  obtained  easily  from  I(Ai,A2) 
by  the  Jacobian  transformation.  The  Fisher  Information  matrix  1(9i,92)  is  1{9i,92)  = 
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((Jij(0i,  02)))  for  i,  j  =  1  and  2,  where 


Iu(0ud2)  = 
I\2{9\,92)  — 
h2(9i,92)  = 


02(?7i0i  +  T162) 

Wl+02)2 


•^21  (^li  ^2)  — 


(n  —  m) 


(0!+02)2’ 
0i(m02  +  nd  1) 

^PT+^j2"' 


Therefore  if  6  —  (0i,  02)  and  9  =  (0i,  0 2),  then  from  the  asymptotic  theory  of  the  MLE’s,  see 
Miller  (1981),  we  have 

(0-0)-»/V2(  0,I-‘(«i,«2)) 
where  I~!  ($i  >  $2)  =  /“1($l  >  ^2)  for  i,j  =  1  and  2  and 


InM) 


In'(8  i,«j) 

/£'(»!,  *j) 


0f  (7x162  ■+■  n01) 

mn92 

hi(9i,92)  = 

02(m0i  +  ^2) 
mn9\ 


9\92{ji  —  m) 


mn 


Note  that  here  1(0;,  92)  is  the  Fisher  Information  matrix  for  the  whole  sample.  To  obtain  the 
confidence  interval  on  9\,  we  substitute  the  true  value  of  the  parameters  by  the  corresponding 
MLE’s  in  the  expression  of  1(0!,  02). 


4.3:  Bootstrap  Confidence  Intervals: 

In  this  section  we  propose  two  bootstrap  confidence  intervals,  the  percentile  bootstrap 
confidence  intervals  suggested  by  Efron  (1982)  and  the  bootstrap-t  confidence  intervals  sug¬ 
gested  by  Hall  (1988).  Hall  (1988)  showed  that  the  bootstrap-t  confidence  interval  is  better 
than  the  percentile  bootstrap  confidence  intervals  from  an  asymptotic  point  of  view.  Al¬ 
though  the  finite  sample  properties  is  not  yet  known. 

We  propose  the  following  percentile  confidence  interval: 


[1]  From  (xi,  £1),  . . .,  (xm,  £m)  obtain  the  bootstrap  sample  (xj,  £J), . . (x^,  5^)  by  resam¬ 
pling  with  replacement  and  from  the  sample  (xm+i,  *), . . .,  (x„,  *)  obtain  the  bootstrap 
sample  (x^+1,  *),  . . .,  (x*,  *)  again  by  resampling  with  replacement. 

[2]  From  the  bootstrap  sample  (xj,  £*)...,  (x^,<5^,),  (xm+i,*),  . . (xn,*)  obtain  the  boot¬ 
strap  estimates  of  all  the  unknown  parameters.  For  any  unknown  parameter,  say  0, 
denote  the  bootstrap  estimates  of  0  as  0*. 

[3]  Repeat  the  process  1-2  NBOOT  times. 
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[4]  Let  CDF(t)  =  P,(9*  <  t)  be  the  cumulative  distribution  of  9*,  the  bootstrap  estimator 
of  the  parameter  9,  then  from  the  NBOOT  9*  obtain  the  upper  bound  and  the  lower 
bound  of  the  100(1-  2a)  %  bootstrap  confidence  bound  for  9  as  follows.  For  a  given  a 

define  9b00t(a)  =  CDF  (a),  then  the  approximate  100(1-  2a)  %  confidence  interval 
for  9  is  given  by 

0boot  (o:) ,  9f,00t  ( 1  a)) 

We  propose  the  following  algorithm  for  computing  the  bootstrap-t  confidence  intervals: 

[1]  From  . . .,  ( xm ,  6m )  obtain  the  bootstrap  sample  (x*,  <5J), . . (xj^,  5^)  by  resam¬ 

pling  with  replacement  and  from  the  sample  (xm+i,  *), . . ( xn ,  *)  obtain  the  bootstrap 
sample  (xj^+1,  *),  . . (x*,  *)  again  by  resampling  with  replacement. 

[2]  From  the  bootstrap  sample  (xj,  5*). . .,  (xj^,  5^),  (xm+i,  *),  . . .,  (xn,  *)  obtain  the  boot¬ 
strap  estimates  of  all  the  unknown  parameters.  For  any  unknown  parameter,  say  9, 
denote  the  bootstrap  estimates  of  9  as  9*. 

[3]  For  any  unknown  parameter  9 ,  compute 

jm  =  sMt  Z  j) 

a* 

where  9  is  the  MLE  of  9  and  a*  is  the  estimated  standard  error  of  9*. 

[4]  Repeat  the  process  1-3  NBOOT  times. 

[5]  From  the  NBOOT  T*  obtain  the  upper  bound  and  the  lower  bound  of  the  100(1-  2a) 
%  bootstrap-t  confidence  bound  for  9  as  follows.  Let  CDFN( t)  =  P,(9*  <  t )  be  the 
cumulative  distribution  of  T*,  then  for  a  given  a  define 

9boot-t{a)  =  9  +  n~?aCDFN  \a) 

The  approximate  100(1-  2a)  %  confidence  interval  for  9  is  given  by 

( ^boot—t  (a) ,  ^6oo(— t  ( 1  o;) ) 


5.  WEIBULL  FAILURE  DISTRIBUTIONS: 
5.1:  Estimation  of  the  Parameters: 
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In  this  section  we  assume  that  Xji's  are  Weibull  random  variable  with  parameters  (a,  Xj) 
for  j  =  1  and  2  and  for  i  =  1,  2,  . . n.  The  distribution  function  Fj(.)  of  Xji  has  the 
following  form: 


Fj(t)  =  (1  -  e-Ajt°). 


(5.1) 


We  assume  that  the  lifetime  distribution  of  the  different  causes  follow  the  Weibull  distribu¬ 
tions  which  has  different  scale  parameters  but  they  have  the  same  shape  parameter,  which 
is  a  quite  practical  assumption,  see  for  example  Rao  et  al.  (1991).  We  introduce  one  more 
parameter  in  the  model,  which  gives  more  flexibility  in  the  hazard  rate,  which  was  constant 
in  case  of  the  exponential  distribution.  The  hazard  rate  or  the  instanteneous  dath  rate  due 
to  cause  1  is  given  by 


gUQ 

m 


aXrf*'-1, 


and  the  mean  lifetime  due  to  cause  1  is 


E(xu)  =  A-r  (i  +  i) . 

Af  V 

It  is  well  known  that  it  can  be  increasing  or  decreasing  depending  on  whether  a  >  1  or  a.  <  1. 
For  a  =  1,  it  has  a  constant  hazard  rate,  because  it  becomes  an  exponential  distribution. 
The  relative  risk  rate,  7r,  due  to  cause  1  is 


r  OO 

7T  =  P(XU  <  X2i)  =  /  aiX  1*°"1e"A,sVAa*a<to 
J  0 


Al 

At  +  A2 


which  is  independent  of  the  shape  parameter  a  and  it  is  the  same  as  the  exponential  case. 
The  log  likelihood  function  of  the  observed  data  as  given  in  Section  3,  becomes 


ln(L)  =  nln(a)  +  r\ln(Xi)  +  r2/n(A2)  +  (a  -  1)  ln(Xi)  -  (Ai  +  A2)  ^2xf 
+  (n  —  m)ln(  At  +  A2). 


i= 1 


t—  1 


Then  taking  the  derivatives  with  respect  to  the  unknown  parameters  a,  Ai  and  A2  and 
equating  them  to  zeros,  we  get 


Ai(a)  = 


n 


n 


A2(a)  = 


n 


m  £”=1  x? 


We  put  Ai  and  A2  in  the  expression  of  In  (L)  above  and  maximize  with  respect  to  a.  We  do 
not  have  any  explicit  expresion  for  d.  We  obtain  d  by  maximizing 


ln[L(pt )] 


n 

nln(a)  +  riln(Xi(a))  +  r2Zn(A2(o!))  +  (a  -  1)  M^i) 

i=l 

(Ai(a)  +  A2(a))  +  (n~  m)ln(\i(a)  +  A2(a)), 

:=1 
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with  respect  to  a.  Once  we  obtain  d ,  we  obtain  the  maximum  likelihood  estimators  of  Ai 
and  A2  as  Ai(d)  and  A 2(d)  respectively.  From  the  invariance  principle  of  the  MLE’s  we  can 
say  that  the  MLE’s  of  the  relative  risk  rate  due  to  cause  1  is 


Ai(dQ 

Ai(d)  +  A2(d) 


and  also  the  MLE  of  the  mean  lifetime  due  to  cause  1  is 


n 


For  known  a  the  distributions  of  \i{a)  or  A2(a)  can  be  obtained  similarly  as  Section  3.  But 
if  a  is  unknown  then  the  exact  distribution  of  Ai(d)  or  A2(d)  is  not  possible  to  obtain.  So 
we  have  to  rely  on  the  asymptotic  distribution  only. 


5.2:  Confidence  Intervals: 


In  this  section  we  provide  the  confidence  intervals  of  the  different  parameters.  Since  the 
exact  confidence  intervals  are  not  possible  to  obtain  when  the  shape  parameter  is  unknown, 
we  propose  the  asymptotic  confidence  intervals  and  also  two  different  bootstrap  confidence 
intervals.  The  asymptotic  result  can  be  stated  as  follows: 

(d  —  a,  Ai  —  Ai,  A2  —  A2)  -4  A^(0, 1  l(a:,  Ai,  A2)). 


Here  I(a,  Ai,  A2)  is  the  Fisher  Information  matrix  for  the  parameters  (o,  Ai,  A2).  The  matrix 
I  =  ((Iij))  for  i,j  =  1,2  and  3  are  as  follows: 


/n(o:,  Aj ,  A2) 

7i2(q:,  Ai,  A2) 
hsi®,  Ai,  A2) 

722(q:,  Ai,  A2) 
hai&i  Ai,  A2) 
723(q:,  Ai,  A2) 


n  I  — —  +  (Aj  +  A2)K 
iaz 

nU  =  72i(q:,  Ai,  A2), 
nU  =  731(q:,  A  i ,  A2 ) , 

7tiA2  +  nAj 

Ai  (A  i  +  A2)2  ’ 
mAi  +  nA2 
A2(Ai  +  A2)2  ’ 

(n  —  m) 


(Aj  +  A2)2 


=  ^32(0!,  Ai,  A2). 


Here  U  =  E(Xa  ln(X))  and  V  =  E(Xaln(X).ln(X)),  where  X  is  distributed  as  Weibull  (a, 
(Ai  +  A2)) .  We  propose  to  use  bootstrap  confidence  intervals  similarly  as  subsection  4.3. 
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6.  NUMERICAL  EXPERIMENTS: 


In  this  section  we  present  some  numerical  results  to  see  how  the  different  methods  behave 
for  small  sample  sizes  and  also  for  different  parametric  values.  All  these  numerical  works 
are  performed  at  the  University  of  New  Brunswick  using  the  HP  workstation.  We  use  the 
method  of  Press  et  al.  (1986)  for  the  random  deviate  generator.  We  consider  the  cases  when 
the  lifetimes  are  exponential  and  also  when  the  lifetimes  are  Weibull.  We  mainly  observe 
the  behavior  of  the  MLE’s  in  terms  of  their  biases  and  in  terms  of  their  variances.  We 
also  compare  the  performances  of  the  different  proposed  confidence  intervals  in  terms  of  the 
coverage  percentages  and  also  in  terms  of  their  average  confidence  lengths. 

6.1:  Exponential  Case: 

In  this  subsection  we  present  the  results,  when  the  lifetimes  are  exponential.  Since  for 
the  exponential  lifetime,  the  biases  and  the  asymptotic  variances  are  functions  of  ^  only, 
we  keep  9\  —  1  fixed,  and  consider  different  values  of  62  =  1.25,  1.50,  ...,  2.50.  We  take 
sample  sizes  n  =  10,  20,  30,  40  and  50.  In  all  the  cases  we  assume  10%  of  the  data  are 
incomplete.  We  draw  random  samples  for  different  values  of  n  and  62  and  compute  the 
MLE’s  of  and  62,  we  also  compute  the  four  different  95%  confidence  intervals,  namely 
(1)  asymptotic  (Asymp)  (2)  approximate  (Approx)  (3)  percentile  Bootstar  (Boot-p)  and 
(4)  Bootstrap-t  (Boot-t)  confidence  intervals.  We  replicate  the  process  one  thousand  times 
and  compute  the  average  values  of  the  MLE’s,  the  variances,  the  biases  and  the  absolute 
biases.  For  the  different  confidence  intervals  we  compute  the  coverage  percentages  and  also 
the  average  confidence  lengths.  The  results  of  9\  and  92  are  quite  similar  in  nature  so  we 
present  the  results  only  of  9\.  The  average  values  of  9±,  the  variances,  the  absolute  biases 
and  the  negative  biases  are  reported  in  Table  6.1.  The  average  error  bounds  (half  of  the 
confidence  length)  and  the  corresponding  coverage  percentages  are  reported  within  brackets 
for  different  methods  in  Table  6.2. 

Some  of  the  points  are  very  clear  from  these  experiments.  From  Table  6.1,  it  is  clear  that 
as  sample  size  increases,  the  biases  and  the  variances  decrease.  It  implies  that  the  MLE’s  are 
asymptotically  unbiased  and  they  are  consistent  estimators  of  the  corresponding  parameters. 
From  the  Table  6.1  and  Table  3.1,  it  is  observed  that  the  theoretical  biases  and  simulated 
biases  are  quite  close  to  each  other.  As  92  increases,  equivalently  as  increases,  the  biases, 

the  variances  of  9 1  decrease  and  the  corresponding  biases  and  variances  of  92  increase  (not 
reported  here).  It  is  not  very  surprising,  because  as  increases  the  mean  life  due  to  cause 
1  decreases  compared  to  the  mean  life  due  to  cause  2.  It  is  expected  that  as  increases, 
the  sample  consists  more  deaths  due  to  cause  1  than  due  to  cause  2.  Therefore  the  sample 
has  more  information  about  9i  than  92  . 

From  Table  6.2  it  is  clear  that  as  sample  size  increases  or  increases,  the  average 
confidence  lengths  of  9\  decrease  for  all  the  four  methods.  In  case  of  92  it  is  observed  (not 
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reported  here)  that  as  sample  size  increases  the  average  confidence  lengths  decrease  but  as 
increases,  the  average  confidence  lengths  increase.  It  is  not  very  surprising  as  it  has  been 
mentioned  previously.  Among  the  different  methods,  it  is  clear  that  all  the  methods  work 
quite  well  if  the  sample  size  is  large  and  all  of  them  are  able  to  keep  the  nominal  coverage 
percentage.  Although  for  small  sample  sizes,  mainly  for  n  =  10  all  the  methods  can  not 
maintain  the  nominal  coverage  percentage.  It  is  observed  that  for  ail  the  methods  except 
the  approximate  one,  the  coverage  perentages  are  slightly  lower  than  the  nominal  level, 
where  as  for  the  approximate  method  the  coverage  percentages  are  slightly  higher  than  the 
nominal  levels.  Between  the  approximate  and  the  asymptotic  method,  the  average  lengths  of 
the  asymptotic  methods  are  slightly  lower  than  that  of  the  approximate  method.  Comapring 
the  percentile  bootstrap  method  and  the  bootstrap-t  method,  it  is  observed  that  bootstrap- 
t  is  preferred  in  terms  of  the  confidence  lengths  although  their  coverage  percentages  are 
almost  equal  in  all  the  cases  considered.  Both  the  methods  can’t  maintain  the  coverage 
percentages  at  least  for  small  samples.  Now  comparing  the  computational  complexities,  the 
asymptotic  method  is  the  easiest  to  obtain.  Approximate  confidence  interval  can  be  obtained 
by  equating  two  non-linear  equation  where  as  the  bootstrap  methods  can  be  obtained  by 
resampling  from  the  original  sample.  It  is  observed  that  bootstrap  methods  take  longer 
times  than  the  approximate  method.  Considering  all  the  points  it  is  recommended  that  for 
small  sample  sizes,  approximate  confidence  interval  can  be  used  where  as  for  large  sample 
the  asymptotic  confidence  bound  is  preferred. 

6.2:  Weibull  Case: 

For  the  Weibull  lifetime  distribution,  we  mainly  consider  the  MLE’s  of  the  c*’s  and  A’s  for 
different  sample  sizes.  We  consider  a.  =  1,  Ai  =  1,  A2  1  =  1.25,  1.50,  1.75,  2.00,  2.25  and  2.50. 
We  take  n  =  10,  20,  30,  40  and  50.  For  a  particular  choice  of  a,  n  and  A2  we  draw  a  random 
sample  from  Weibull  lifetime  distribution  and  compute  the  MLE’s  of  Aj,  A2  and  a  when  10% 
data  are  incomplete.  We  also  compute  the  three  different  95%  confidence  bounds  namely 
the  asymptotic  one,  the  percentile  bootstrap  and  the  bootstrap-t  confidence  intervals.  We 
replicate  the  process  one  thousand  times  and  compute  the  average  biases,  absolute  biases 
and  the  variances  of  the  MLE’s  over  one  thousand  replications.  The  average  variances,  the 
negative  bises  and  the  absolute  biases  of  Ai  are  reported  in  Table  6.3.  The  average  error 
bounds  and  the  coverage  percentages  (within  brackets)  are  reported  in  Table  6.4.  Since  the 
results  of  A2  and  a  are  quite  similar  to  that  of  Ai  they  are  not  reported  here. 

Some  of  the  points  are  very  clear  from  this  experiment.  From  the  Table  6.3  it  is  clear  that 
as  sample  size  increases  the  variances  the  biases  and  the  absolute  biases  all  decrease  which 
indicates  the  consistency  of  the  MLE’s  even  in  the  Weibull  case  also.  Comparing  Table  6.1 
and  Table  6.3  it  is  observed  that  the  biases,  variances  and  the  absolute  biases  of  9\  are  less 
than  the  corresponding  biases,  variances  and  the  absolute  biases  of  Ai,  although  both  are 
estimating  the  same  parameter  Ai  =  flf1  =  1  in  this  case.  Which  is  not  very  surprising 
because  in  the  exponential  case  only  two  parameters  to  estimate  whereas  for  the  Weibull 
case  there  are  three  unknown  parameters. 
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From  Table  6.4,  it  is  clear  that  none  of  the  methods  are  able  to  maintain  the  coverage 
percentages  for  small  sample  sizes.  Although  for  the  large  sample  sizes  all  the  three  methods 
behave  reasonably  well.  As  far  as  the  confidence  lengths  are  concerned  the  confidence  lengths 
due  to  the  asymptotic  methods  have  the  marginally  smallest  size  than  the  other  two  and 
coverage  percentages  are  close  to  the  nominal  value  at  least  for  large  sample  sizes.  Therefore 
for  moderate  or  large  sample  sizes  asymptotic  method  can  be  used. 


7.  DATA  ANALYSIS: 


In  this  section  we  consider  one  real  life  data  set  from  Lawless  (1982,  page  491).  It  consists 
of  failure  or  censoring  times  for  36  appliances  subjected  to  an  automatic  life  test.  Failures 
were  classified  into  18  different  modes,  though  among  33  observed  failures  only  7  modes  are 
present  and  only  modes  6  and  9  appear  more  than  once.  We  are  mainly  interested  in  the 
failure  mode  9.  The  data  consist  of  two  causes  of  failure  (5  =  1  (failure  mode  9)  and  6  =  2 
(all  other  failure  modes.  The  data  are  given  below,  which  indicates  the  failure  times  and  the 
cause  of  failure  if  available. 

Data  Set:  (11,2),  (35,2),  (49,2),  (170,2),  (329,2),  (381,2),  (708,2),  (958,2),  (1062,2),  (1167,1), 
(1594,2),  (1925,1),  (1990,1),  (2223,1),  (2327,2),  (2400,1),  (2451,2),  (2471,1),  (2551,1),  (2565,*), 

(2568.1) ,  (2694,1),  (2702,2),  (2761,2),  (2831,2),  (3034,1),  (3059,2),  (3112,1),  (3214,1),  (3478,1), 

(3504.1) ,  (4329,1),  (6367,*),  (6976,1),  (7846,1),  (13403,*) 

Here  we  have  n  =  36,  m  =  33,  r i  =  17,  r2  =  16,  =  99245.  Therefore  using  the 

exponential  lifetime  distribution,  the  ML  estimators  9\  =  5351.45  and  62  =  5685.91.  The 
estimates  of  the  mean  life  due  to  cause  1  and  cause  2  become  5351.45  and  5685.91  respectively. 
The  ML  estimators  of  the  relative  risk  rate  due  to  cause  1  is  7r  =  .5152  and  due  to  cause  2 
is  1  -  7T  =  .4848.  The  following  95  %  confidence  intervals  are  obtained  for  0\  and  62  by  using 
different  methods. 


8 

1 

8 

2 

Methods 

LB 

UB 

LB 

UB 

Asymp. 

2862.73 

7840.16 

2956.68 

8415.14 

Approx. 

4451.45 

6251.28 

4785.91 

6585.22 

Boot-p 

3683.07 

7407.33 

2464.83 

9406.65 

Boot-t 

3534.64 

6992.28 

3421.01 

7103.94 

Using  Weibull  lifetime,  we  obtain  the  estimates  of  6\  =  1  =  6980.22,  §2  =  A2 1  = 

7416.49  and  &  =  1.0321.  The  95  %  confidence  band  of  a  becomes  (.7625,  1.3016).  Since  this 
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interval  contains  one,  therefore  we  cannot  reject  the  null  hypothesis  Hq:  a  =  1  at  the  95  % 
level.  Therefore  we  conclude  that  the  competing  lifetimes  are  exponential  and  they  have  the 
constant  failure  rates.  Since  the  simulation  results  indicate  that  the  approximate  method 
works  quite  well  for  the  exponential  case,  we  use  the  approximate  confidence  bands  for  the 
parameters  6 1  and  d2  respectively,  which  are  given  in  Table  7.1. 


8.  CONCLUSIONS: 


In  this  paper  we  consider  estimation  of  the  parameters  of  the  competing  risks  model 
when  the  data  may  not  be  complete.  We  consider  two  different  lifetime  distributions  of  the 
competing  causes,  namely  exponential  and  Weibull.  We  obtain  the  exact  distribution  of  the 
MLE’s  of  the  mean  lifetime,  when  the  lifetime  distributions  are  assumed  to  be  exponential. 
We  propose  approximate  confidence  bands  for  the  mean  lifetime  and  compare  their  perfor¬ 
mances  with  the  asymptotic  confidence  bands  and  two  other  bootstrap  confidence  bands. 
It  is  observed  that  approximate  confidence  bands  work  quite  well  for  the  exponential  case. 
When  the  lifetime  distributions  are  Weibull,  it  is  observed  that  MLE’s  behave  reasonably  well 
and  similarly  as  the  exponential  case  they  also  provide  consistent  estimates  of  the  unknown 
parameters.  To  obtain  the  confidence  bounds  of  the  unknown  parameters  the  asymptotic 
results  can  be  used  for  moderate  or  large  sample  sizes  but  for  small  sample  sizes  more  work 
is  needed. 

Another  important  aspect  which  is  not  addressed  here  is  the  analysis  when  the  competing 
risks  may  not  be  independent  and  when  some  of  the  causes  of  failure  are  not  known.  It  is 
a  difficult  problem.  One  way  it  can  be  handled  through  mixture  model  formulation  as  was 
suggested  by  Babu  et  al.  (1992)  or  Kundu  et  al.  (1992).  It  will  be  reported  else  where. 
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Table  6.1 

Average  values  of  9\,  Variance  of  9\ ,  bias  and  absolute  bias  of  the  LSE’s  when  10%  of  the 
date  are  incomplete.  The  lifetimes  are  exponential. 


n 

Average 

e2  =  1.25 

02  =  1-50 

02  =  2.25 

02  =  2.50 

10 

«, 

V(«,) 

Bias 

|Bias| 

1.1261 

.4107 

.1261 

.4107 

1.1018' 

.3459 

.1018 

.3801 

1.0947 

.2872 

.0947 

.3575 

1.0741 

.2409 

.0741 

.3343 

1.0776 

.2227 

.0776 

.3235 

1.0615 

.1902 

.0615 

.3056 

20 

V(«.) 

Bias 

|Bias| 

1.0411 

.1410 

.0411 

.2510 

1 

1.0303 

.0894 

.0303 

.2273 

1.0291 

.0850 

.0291 

.2223 

30 

fx 

vA) 

Bias 

|Bias| 

1.0243 

.0636 

.0243 

.1953 

1.0228 

.0581 

.0228 

.1854 

1.0173 

.0523 

.0173 

.1782 

40 

V(«.) 

Bias 

|Bias| 

1.0264 

.0564 

.0264 

.1816 

1.0215 

.0500 

.0215 

.1727 

1.0134 

.0440 

.0134 

.1647 

1.0175 

.0437 

.0175 

.1625 

1.0095 

.0402 

.0095 

.1558 

1.0097 

.0373 

.0097 

.1522 

50 

vA) 

Bias 

|Bias| 

1.0192 

.0425 

.0192 

.1588 

1.0151 

.0388 

.0151 

.1534 

1.0153 

.0362 

.0153 

.1483 

1.0116 

.0339 

.0116 

.1438 

1.0078 

.0298 

.0078 

.1360 

20 


Table  6.2 

Different  confidence  intervals  of  9\,  when  10%  of  the  date  are  incomplete.  The  lifetimes  are 
exponential.  The  average  length  of  the  error  bound  and  the  corresponding  coverage 
probability  (within  bracket)  are  reported.  The  nominal  level  is  95%. 


n 

Methods 

02  =  1.25 

02  =  1.50 

02  =  1.75 

02  =  2.00 

02  =  2.25 

Asymp 

1.089(.917) 

.979(.900) 

.949(.898) 

.908(.897) 

.892(.867) 

.865(.847) 

10 

Approx 

.703(.928) 

.707(.938) 

.701(.925) 

.697  (  .931) 

.692  (.909) 

Boot-p 

1.528(.901) 

1.442(.894) 

1.298(.863) 

1.219(.870) 

1.166(.873) 

Boot-t 

1.140(.898) 

1.076(.891) 

.985(.873) 

.948(.867) 

.888(.862) 

.919(.868) 

Asymp 

.656(.917) 

.627(.923) 

.596(.908) 

.576(.911) 

.559(.940j 

.552(.918) 

20 

Approx 

.625(.952) 

,601(.962) 

.605(.940) 

.590  ( .962) 

.522(.956) 

.523(.955) 

Boot-p 

.969(.933) 

.841(.939) 

.774(.923) 

.683  (  .935) 

.691(.931) 

.668(.922) 

Boot-t 

.794(.935) 

.719(.931) 

.683(.926) 

,617(.925) 

.622(.928) 

.608(.924) 

Asymp 

.516(.913) 

.492(.918) 

.471  (.928) 

.462(.933) 

.454(.934) 

.441(.938) 

30 

Approx 

.547(.952) 

.529(.952) 

.508(.960) 

.494  (  .956) 

.485(.948) 

Boot-p 

.645(.940) 

.604(.951) 

.561(.944) 

.535  (  .932) 

.514(.944) 

lai 

Boot-t 

.585(.939) 

.558(.951) 

.528(.934) 

.506(.928) 

.490(.945) 

Asymp 

.444(.928) 

.427(.944) 

.408(.924) 

.400(.946) 

.392(.940) 

.377(.928) 

40 

Approx 

.473(.950) 

.449(.964) 

.430(.953) 

.424(.963) 

.417(.953) 

.408(.956) 

Boot-p 

.524(.944) 

.503(.938) 

.466(.945) 

.452  ( .945) 

.435(.950) 

.431  (.944) 

Boot-t 

.492(.935) 

.477(.941) 

.447(.941) 

.436(.939) 

.421(.944) 

.419(.942) 

50 

Asymp 

.399(.946) 

,378(.931) 

.362(.948) 

.352(.933) 

.345(.931) 

.340(.939) 

Approx 

.424(.966) 

.399(.963) 

.384(.949) 

.372  (.948) 

.365(.957) 

.359(.967) 

Boot-p 

.466(.962) 

.438(.948) 

.415(.960) 

.401  (  .953) 

.393(.949) 

.386 (.948) 

Boot-t 

.445(.965) 

.421(.947) 

.404(.961) 

.389  (  .955) 

.383(.949) 

.369 (.948) 

21 


Table  6.3 

Average  values  of  Ai,  Variance  of  A1}  bias  and  absolute  bias  of  the  LSE’s  when  10%  of  the 

date  are  incomplete.  The  lifetimes  are  Weibull. 


n 

Average 

A2 1  =  1.25 

AJ 1  =  1.50 

Aj 1  =  1.75 

AJ1  =  2.25 

10 

Al 

V(A.) 

Bias 

|Bias| 

1.2356 

.8067 

.2356 

.5006 

1.2138 

.5599 

.2138 

.4530 

1.1926 

.6956 

.4926 

.4334 

1.1994 

.5110 

.1994 

.4315 

1.1925 

.8079 

.1925 

.4260 

1.1385 

.3775 

.1385 

.3742 

20 

A, 

V(A,) 

Bias 

|Bias| 

1.0996 

.1553 

.0996 

.2859 

1.1116 

.1318 

.1116 

.2769 

1.0749 

.1250 

.0794 

.2696 

1.0844 

.1278 

.0844 

.2649 

1.0972 

.1077 

.0972 

.2499 

1.0781 

.1241 

.0781 

.2556 

30 

Ai 

V(A.) 

Bias 

|Bias| 

1.0499 

.0793 

.0499 

.2117 

1.0426 

.0713 

.0426 

.2030 

1.0448 

.0623 

.0448 

.1962 

1.0488 

.0661 

.0488 

.1982 

1.0591 

.0764 

.0591 

.2059 

1.0357 

.0627 

.0357 

.1925 

40 

Ai 

V(Ai) 

Bias 

|Bias| 

1.0425 

.0542 

.0425 

.1819 

1.0429 

.0502 

.0429 

.1755 

1.0452 

.0470 

.0452 

.1670 

1 

1 

1.0293 

.0457 

.0293 

.1672 

1.0250 

.0410 

.0250 

.1591 

50 

Ai 

V(A!) 

Bias 

|Bias| 

1.0199 

.0404 

.0199 

.1595 

1.0274 

.0408 

.0274 

.1577 

1 

1 

1.0267 

.0372 

.0267 

.1502 

I 

m  u 

1  1 

22 


Table  6.4 

Different  confidence  intervals  of  A^  when  10%  of  the  date  are  incomplete.  The  length  of 
the  confidence  intervals  and  the  corresponding  coverage  probability  (within  bracket)  are 
reported.  The  lifetime  distributions  are  Weibull. 


n 

Methods 

Aj  1  =  1.25 

AJ1  =  1;50 

II 

-<i 

cn 

Aj  1  =  2.00 

AJ1  =  2.25 

10 

Asymp 

Boot-p 

Boot-t 

1.131(.929) 

1.282(.875) 

1.152(.865) 

1.063(.932) 

1.273(.854) 

1.121(.840) 

1.012(.925) 

1.176(.848) 

1.052(.834) 

.993  (.931) 
1.136(.836) 
1.047(.832) 

.973(.893) 

1.087(.841) 

.972(.831) 

.893(.890) 

1.026(.838) 

.892(.806) 

Asymp 

Boot-p 

Boot-t 

.610(.932) 
.896  (.919) 
.699  (.911) 

.593(.948) 
.824  (.923) 
.657(.921) 

.575(.939) 
.802  (.925) 
.609  (.913) 

30 

Asymp 

Boot-p 

Boot-t 

.521(.938) 
.694  (.935) 
.602  (.938) 

.496(.941) 
.634  (.951) 
.571  (.943) 

.479(.956) 
.605(.945) 
.544 (.938) 

.471(.953) 
.598  (.931) 
.540  (.943) 

.467(.937) 

.572(.950) 

.529(.944) 

.449(.942) 

.569(.939) 

.514(.941) 

Asymp 

Boot-p 

Boot-t 

.446(.951) 
.554  (.934) 
.516  (.929) 

.427(.952) 
.527  (.940) 
.485  (.942) 

.416(.952) 
.497  (.942) 
.474  (.937) 

.400(.948) 
.493  (.941) 
.455  (.943) 

.392(.947) 

.476(.947) 

.453(.939) 

.383(.943) 

.462(.940) 

.437(.945) 

50 

Asymp 

Boot-p 

Boot-t 

.392(.957) 
.470  (.963) 
.454  (.942) 

.377(.952) 
.453  (.944) 
.429  (.940) 

.365(.943) 
.428  (.959) 
.417  (.939) 

.356(.941) 
.423  (.952) 
.416  (.946) 

.348(.950) 

.407(.948) 

.399(.947) 

.341(.959) 
.406(.950) 
.391  (.948) 

23 


